The thermopower as a signature of quantum criticality in heavy fermions 
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We present a series of arguments showing that the Seebeck coefficient can be used as a deci- 
sive experiment to characterize the nature of the quantum critical point (QCP) in heavy fermion 
compounds. Being reactive almost exclusively to the presence of delocalized entropic carriers, the 
Seebeck coefficient shows a drastic collapse at the Kondo breakdown QCP, as the reconstruction of 
the Fermi surface takes place. In contrast, around a Spin Density Wave (SDW) QCP, the Seebeck 
coefficient is broadly symmetric. We discuss the possibility of a change of sign at the QCP, the 
characteristic variation of \S/T\ with temperature and external parameter, as well as the capacity 
of the Seebeck coefficient to distinguish between localized and itinerant anti-ferromagnetism (AF). 
Suggestions of experiments are given in the case of four non conventional compounds : YbRli2Si2, 
Ce(Mn)In 5 , CeCiie-zAu^ and URu 2 Si 2 . 



INTRODUCTION 

In the mid-eighties, heavy fermion compounds were in- 
tensively studied for their heavy Fermi liquid properties 
[IHH . Although those dense rare earth lattices were made 
of big localized atoms supporting big magnetic moments, 
the low energy properties remained in the universality 
class of the Landau Fermi liquid theory of metals, with a 
characteristic low temperature saturation of the Sommer- 
feld coefficient 7 = C/T of the magnetic susceptibility 
Xo ~ Gst and a Independence of the corrections to the 
residual resistivity p(T) — po = AT 2 . The effective mass 
of the Landau quasi-particles is strongly re-normalized, 
up to a factor of ~ 1000 for UBei3, but still the obser- 
ved low temperature properties didn't depart from the 
Landau theory of metals. 

The situation changed drastically in the mid-nineties, 
with the observation that, under the application of an 
external tuning parameter like chemical doping, pressure 
or magnetic field, the specific heat coefficient doesn't sa- 
turate when the temperature is lowered Q • This anoma- 
lous property was rapidly attributed to the presence of a 
QCP where the system orders anti-ferromagnetically at 
(theoretically) vanishing temperature. The strong quan- 
tum fluctuations induced at the vicinity of a zero tem- 
perature phase transition were suggested to be respon- 
sible for this violation of the Landau theory of metals. 
Rapidly other properties were shown not to follow the 
universal Landau paradigm. In many compounds the re- 
sistivity is linear, or quasi-linear in temperature over two 
or three decades in energy [l|, [5| , the magnetic susceptibi- 
lity shows some anomalous exponents with temperature 
like in CeCus.gAuo.iQ. As the tuning parameter evolves 



from the heavy Fermi liquid phase towards the QCP, the 
A-coefficient of the T 2 -resistivity shows a divergent trend 
with respect to the tuning parameter [7]. The effective 
band mass, shown by de Haas van Alphen experiments 
is strongly re- normalized at the approach of a QCP0, 
which is a remarkable fact since the renormalization of 
the band mass is essentially due to elastic scattering pro- 
cesses. All over the years, an average of 20 compounds 
were found to have anomalous physical properties, when 
fine tuning with an external parameter was performed. 
Those findings are well summarized in various review ar- 
ticles [2j, y, l9( and we refer the reader to them for further 
details. 

In heavy fermion compounds the study of the thermo- 
power started three decades ago [n| and a few syste- 
matic features were already clarified. The high tempera- 
ture thermopower is typically large due to the interplay 
of incoherent spin fluctuations and crystal field effects. 
Like most of the thermodynamic and transport proper- 
ties, the termopower shows a maximum corresponding 
to the bandwidth of the f-electrons. This scale is some- 
times refered to as the lattice coherence scale To [ll| . 
The sign of the Seebeck coefficient has been shown to 
depend crucially on the position of the f-resonance level 
with respect to the Fermi surface. In the Ce (4f 1 )-series 
the f-resonance sits above the Fermi level, which leads to 
a positive Seebeck coefficient, while in the case of the Yb 
(4f 13 )-based compounds, the f-level is below the Fermi 
level which leads to a negative Seebeck coeffcicient. The 
case of the U (4f 2 )-based compounds is more controver- 
sial, since some compounds are compensate metals with 
very low carrier density, like URu2Si2 and have a nega- 
tive Seebeck coefficient (12| while for UPt3 no detectable 
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signal has been observed [l3j]. For UBei3 , S/T is strong 
and negative Q while for UPd 2 Al 3 and UNi 2 Al 3 [l]| it 
is small and positive, revealing the complexity of the f- 
electron structure in U-based compounds. 

Quite remarkably, at low temperature in heavy fer- 
mions, the Seebeck coeffcient divided by the temeprature 
S/T was shown to form a constant ratio with the Som- 
merferld specific heat coefficient 7 [l(| 
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where Nav is the Avogadro number and e the electron's 
charge. This ratio is close to ±1 for the majority of heavy 
fermion compounds. This quasi-universal behavior was 
explained by the observation that although many bands 
are present in a typical band structure of heavy fermions, 
the Seebeck coefficient is mostly sensitive to the position 
of the heaviest band, namely the one with the biggest 
f- character [fT). Since the Sommerfeld ratio is precisely 
sensitive to the heaviest band as well, a quasi-universal 
behavior is to be expected. The formula (fT]) tells us that, 
in the Fermi liquid regime, the thermopower probes the 
specific heat per electron. This ratio can be compared 
with other quasi-universal ratios studied in heavy fer- 
mion systems. The Wilson ration xll l2fl| °f t ne ma ~ 
gnetic susceptibility to the Sommerfeld coefficient and 
the Kadowaki-Wood ratio A/7 2 [2l( of the coefficient A 
of the T 2 resistivity in metals show some universal ratio 
insensitive to the mass renormalization in of the heavy 
fermi liquid. 

In simple metals, the thermoelectric effects are very 
sensitive to the type of scattering involved. In addition 
to the diffusion Seebeck coefficient, the electron-phonon 
interaction produces a phonon-drag component which do- 
minates the behavior in many metals [22|. In the pre- 
sence of various types of scattering the thermopower is 
the sum of the contribution of each scattering process, 
weighted by the resistivity, a rule reminiscent of the Mat- 
thiessen rule for the addition of resistivity, referred to as 
the Nordheim-Gorter rule 



S = 
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In the case of multi-bands systems, the Mott [23] rule 
applies where the Seebeck coefficient for each band is 
weighted by the conductivity 



S = 
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Little is known about the Seebeck coefficient close to 
a QCP. Preliminary studies for CeCue-zAor 24j and 
Ce( Nii_ a; Pd : i;)2Ge2[25| show that the presence of a QCP 
modifies low temperature dependence of the Seebeck co- 
efficient. Two recent studies under magnetic field show 
some striking similarity between thermoelectric effects in 



CeCoLri5 [26| and URu2Si2 (27j . In particular both system 
show a pronounced anisotropy in their thermoelectric res- 
ponse. Lastly, a recent experiment on YbRh 2 Si2 under a 
small magnetic field shows some drastic variations of the 
magnitude of the Seebeck coefficient on both sides of the 

qcp HU. 

Even fewer theoretical studies are available 17 In 
the case of the Spin Density Wave (SDW) QCP, the au- 
thors of Ref.[29| have shown that at the QCP S/T has 
the same variation with temperature as the Sommerfeld 
coefficient 7(T).The low temperature correlation between 
S/T and 7 survives close to a QCP. 

In this discussion paper, we want to address the re- 
levance of thermoelectric properties close to the Kondo 
breakdown QCP. In the next section we give an over- 
view of the Kondo breakdown theory and explain why the 
Seebeck coefficient might be the best probe to characte- 
rize the Kondo breakdown QCP. We make as well some 
distinctions between the SDW scenario and the Kondo 
breakdown, which can lead to experimental discrimina- 
tion between the two QCP. In the next section, we review 
the unconventional properties of QCPs in four heavy fer- 
mion compounds and suggest useful thermoelectric ex- 
periments susceptible to unravel the true nature of the 
QCP. 



UNCONVENTIONAL QCPS : THE KONDO 
BREAKDOWN MODEL 



In this section we review the various QCPs that have 
been suggested to explain the very unconventional be- 
havior observed in heavy fermions. Heavy fermions are 
heavy metals made of big magnetic atoms interacting 
hybridized to a bath of conduction electrons. Many com- 
pounds exhibit magnetic phases, anti-ferromagnetic (AF) 
or frustrated. It was natural to attribute the anomalous 
properties of those compounds to the proximity to a ma- 
gnetic phase transition at T = 0. At this QCP, the Fermi 
surface is destabilized by spin density waves. This scena- 
rio is also called the SDW theory. It has been derived by 
Hertz and revived by Millis [3l| . At the heart of this 
theory is the itinerant character of conduction electrons. 
When a bosonic mode of the type of a SDW interacts 
with conduction electrons, the particle-hole continuum 
produces Landau damping —ico/q, where q is the modu- 
lus of the scattering vector. If the QCP sits at the brink 
of uniform ordering, like for example in the case of a fer- 
romagnet we are in the regime where q — > and \uj\ < q. 

For incommensurate or AF order though, the ordering 
wave vector is finite and the damping takes the form 
—iuj/Q*, where Q* is the modulus of the ordering wave 
vector. In this case the spin susceptibility in the vicinity 
of the QCP writes 



V + r 
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where 7 = 27tto and m is the band mass of the conduc- 
tion electrons. £ is the correlation length which depends 
both on the temperature and on the distance from the 
QCP ; at the QCP £ — > 00. We see that in this model, 
the fluctuations in the imaginary time, also called the 
quantum fluctuations, scale like u> ~ q 2 which defines 
the dynamical exponent z = 2 [32]. The treatment of 
the Hertz-Millis theory requires to integrate the fermions 
out of the partition function, which is an uncontrolled 
operation. A better treatment is given by writing a set of 
self-consistent equations for the polarization and the self- 
energy and using the Migdal theorem to neglect vertices. 
This amounts to performing an Eliashberg treatment of 
this theory 33|, 3J|. Those two techniques give the same 
result, and in the absence of a reliable bosonization of 
the SDW model (35|, they constitute the state of the art. 

The results obtained from the SDW model are sum- 
marized in the left panel in FigTTaland in Table Q] In the 
quantum critical regime, the Sommerfeld coefficient di- 
verges logarithmically at low temperature in dimension 
two but not in dimension three. Its generic scaling with 
the temperature goes like 7 ~ y(«J-2)/2_ rp^g correc tions 
to the electrical resistivity vary like p — po ~ T d / 2 . This 
power law has to be understood as a correction to the 
residual resistivity. It is valid when p — po -C pq. The 
static staggered spin susceptibility varies like T d / 2 . The 
Seebeck coefficient divided by the temperature [29] varies 
like the Sommerfeld coefficient, as T^ 2 ^ 2 . When cros- 
sing the QCP by decreasing or increasing the external 
parameter "x" (here "x" represents pressure, doping or a 
small magnetic field) a doubling of the Brillouin zone is 
observed, but with conservation of the Luttinger theo- 
rem ; the number of electronic carriers is conserved. 

In contrast with the SDW model, the Kondo Break- 
down QCP is not associated with a T = magnetic phase 
transition but with the localization of the f-electrons un- 
der very strong onsite Coulomb potential U. When the 
QCP is crossed from the right to the left by varying "x", 
a complete reconfiguration of the Fermi surface is ob- 
served ; the Fermi surface of the f-electrons becomes hot 
and on the left side of the phase diagram the f-electrons 
don't participate to the transport. This transition has 
been dubbed in some works selective Mott transition, 
the word selective referring to the localization of the 
f-electron while the conduction electrons remain itine- 
rant. The T = phase transition in this case is a de- 
confinement transition for the heavy electron. 

At the QCP the heavy electron splits into three parts, 
(i) the conduction electrons, (ii) the spinons (xa, xV) 
carrying spin and (iii) the holons (b,w) carrying charge, 
describing the break-up of the f-electron at the Mott 
localization. In a field theory language this transition 
is described condensation of the holon-operator 

b = 6 T = bo, within a spontaneous symmetry breaking 
also called Anderson-Higgs transition. Fictitious gauge 
fields are generated to sustain the U(l) gauge symme- 



try. The whole description with spinons and holons can 
be understood as a field-theoretical way of tracking the 
Mott transition, analogous to what was implemented for 
the single band Hubbard model in the early days of the 
cuprate superconductors [36]. 

One of the main differences between the SWD model 
and the Kondo Breakdown is that the Kondo Breakdown 
has a z — 3 quantum critical regime (instead of the z = 2 
quantum critical regime of the AF SDW model) [37| • In 
this regime the typical form of the holon propagator reads 

M 



D-\ q ,u)^^ +q 2 +i 
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where a is a dimensionless number much smaller than 
one, which represents the ratio between the %-spinons 
and the c-electrons bandwidths. As a result, the critical 
exponents are different from the ones of the SDW thoery. 
The Sommerfeld coefficient varies like 7 ~ T^- 3 ^ 3 . The 
resistivity varies like p—po ~ T d l 3 and the Seebeck coeffi- 
cient over the temperature varies like S/T ~ T^ 3 ^ 3 . It 
is worth noticing that in dimension three, the resistivity 
is quasi-linear in temperature, with Ap ~ TLog(T/E*). 
The quasi-linear temperature exponent is not correction 
to the residual resistivity, but a robust exponent due 
to the specificity of the Kondo Breakdown to have two 
kinds of particles, light conduction electrons and almost 
localized spinons. The light conduction electrons scatter 
through the local network of spinons, the scattering pro- 
cess involves the z = 3 critical bosons, producing a trans- 
port exponent quasi-linear in temperature |41| . The cri- 
tical exponents for the SDW and the Kondo Breakdown 
models are compared in Tabl^U The Kondo Breakdown 





SDW 


KBreakdown 


Ap 


rpd/2 


rpd/3 
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ji(d-2)/2 


rp(d-3)/3 


S/T 


jr(d-2)/2 


j>(d-3)/3 



Table I: comparison of the critical exponents for the 
SDW (z = 2) and the Kondo Breakdown (z = 3) 

models, for the resistivity, the Sommerfeld coefficient 
and the ratio of the Seebeck coefficient with the 
temperature. 

model also differs from the SDW theory from because of 
the emergence of an additional scale -called E*-, at the 
QCP.The scale E* , in this model, is due to the presence 
of two types of fermions, the \ spinons and the conduc- 
tion electrons. The two corresponding Fermi surfaces are 
not necessarily close to each other. In the case where they 
are centered, the mismatch q* = \kp — kp\ between the 
Fermi wave vectors produces an additional energy scale 



E* ~ 0.1( 9 7fc^) 3 T 
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Quantum Critical Regime 

z=3 ^ ' ' 

/ 

/ 

/ 

/ heavy Fermi liquid 

QCP x 

Kondo breakdown phase diagram for heavy 
fermions 

FIGURE 1: This figure compares the phase diagram of the Kondo Breakdown QCP with the one of the SDW model 
at and AF phase transition. In the SDW model the f-electrons remain itinerant while crossing the QCP while in the 
Kondo Breakdown, the f-electrons localize below E* on the left-hand side of the phase diagram. In the Kondo 
Breakdown, the scale To marks the onset of a spin liquid, where entropy is quenched by entanglement of the 
f-momcnts with no long range order ; in the SDW To can be any scale associated with the mean-field formation of 

Kondo singlets. A crucial difference between the two QCPs is that the quantum critical regimes have different 
dynamical exponents, with z = 2 for the SDW and z = 3 for the Kondo breakdown. That means that experimental 
observable have different exponents in this regime. The scale E* is typical of the Kondo breakdown, it marks the 
end of the 2 = 3 regime. The regime below E* has Fermi liquid characteristics. AF long range order can occur at the 
vicinity of the Kondo Breakdown QCP but is directly tight with the zero temperature phase transition. It has been 
omitted in this figure. Note that the scale E* made to stop at the boundary of the heavy Fermi liquid regime, since 

this scale is difficult to observe in this regime. 




where To is the scale above which the entropy i?ln2 is 
quenched. In the case where the two Fermi surfaces are 
not centered, the holons condense at finite qo in order to 
re-center them [371). Note that the power 3 exponent in 
Eqn.([n]) ensures that for quite a number of compounds 
E* is small ; typically E* mm < E* < E* max with E* max ~ 
250mK for T ~ 20 K and q*/k F ~ 0.5 and E* mm ~ 
for T ~ 20K and q*/k F ~ OCT 2 . 
The scale E* is a key feature of the Kondo breakdown 
QCP. Below E* the particle-hole continuum is gapped 
out, hence the order parameter reduces to a free boson 
mode below the gap with the dispersion D' b 1 (q,uj) cx 
w + q 2 + £~ 2 - This mode doesn't lead to any appreciable 
contribution to the thermodynamics and transport, and 
the regime below E* can be characterized by small cor- 
rections to the Fermi liquid theory [3{|. The reconfigu- 
ration of the Fermi surface at the QCP can be found in 
Ref. [!(| • The multi-scale character of the QCP, as well as 
the z = 3 regime can be found in the work of Ref.[4l|. 
The Kondo breakdown QCP has already been the object 
some scrutiny by various groups ji^- In particular, two 
DMFT studies are now confirming its existence 4a, |46( . 
The Kondo breakdown can be described through an ef- 
fective low energy field theoretical Lagrangian, which en- 
ables refine the theoretical predictions. The most com- 
plete treatment to date, however, still relies on an Elia- 



shberg theory where the vertices are neglected and the 
self-energies retained. At this stage of development, the 
theory suffers form the fact that the localized spinons are 
described within a fermionic representation of the spin 
(Abrikosov pseudo- fermions). Hence the properties asso- 
ciated with the entropy of the localized spins are poorly 
described. We expect however, that the model gives a 
correct description of the transport properties. 

Another scenario has been proposed in the literature 
to explain the anomalous properties observed in those 
compounds : the locally quantum critical scenario (47l |. 
This theory is also based on a breakdown of the heavy 
Fermi liquid, and thus enter the generic class of " Kondo 
Breakdown" scenarios. However it leads to different re- 
sults as the Kondo breakdown QCP and it is supported 
by a few assumptions that we believe will become experi- 
mentally testable in the near future. The locally quantum 
critical point requires the presence of two dimensional 
spin fluctuations. It predicts some anomalous exponents 
in the spin susceptibility x ~ 7 1-0 - 75 over a wide range 
of the Brillouin zone. Moreover, it is always situated at 
the brink of a magnetic T = phase transition. This pro- 
perty distinguishes it from the Kondo Breakdown QCP, 
which not directly correlated to the occurrence of long 
range magnetic order. 

It is interesting to find out what kind of phase diagram 
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one obtains when the Kondo breakdown and magnetism 
are treated together. At the present moment the avai- 
lable theories don't allow us describe both phenomena 
together in a controlled way, but it is still interesting 
to consider the putative phase diagram one would ob- 
tain. The result is presented in Fig[5J where an additional 
axis of frustration has been added to the system. Due to 
the Rudderman-Kittel-Kasuya-Yoshida (RKKY) interac- 
tions, frustration is naturally generated in the Kondo lat- 
tice, and it is interesting to think that the combined ef- 
fects of crystal fields and geometric frustration vary from 
compound to compound and lead to various magnitude 
and different structure for the AF magnetic order. In the 
3d-phase diagram, a line of Kondo Breakdown QCPs is 
crossing the AF long range order line (at T = 0) at one 
point only, and this crossing is accidental. The line of 
Kondo Breakdown QCPs separates a localized regime on 
the left to an itinerant regime on the right. 

The phase diagram of Fig [2] presents some analogies 
with the two fluids scenario for heavy fermions, which has 
been proposed as a generic mechanism for the formation 
of the heavy Fermi liquid [49]. In order to understand bet- 
ter the analogy, it is instructive to think about those com- 
pounds in terms of entropy. At high temperature, heavy 
fermions are made of a lattice of fluctuating spins, having 
a marginal and incoherent interactions with the conduc- 
tion electrons. As temperature is lowered, this entropy is 
quenched by spin-orrbit and crystal field coupling, but 
still, at intermediate temperature a finite entropy of the 
order of R In 2 per impurity site remains. As temperature 
is lowered further down, two routes open to quench the 
entropy. Either the impurity spins start to entangle toge- 
ther, quenching the entropy via the formation of a para- 
magnetic liquid which is called here a "spin liquid", or the 
entropy can be quenched by the formation of Kondo sin- 
glets, which for the lattice finally leads to the formation 
of the heavy Fermi liquid. Lastly at even lower tempe- 
ratures, long range order of various kinds can achieve 
further quenching of the entropy. Mainly AF order or 
super-conductivity occur as a rule. Those two routes for 
the quenching of the entropy - the formation of a spin 
liquid or the one of a heavy Fermi liquid- have been in- 
tuited forty years ago by Doniach [50||. He was the first 
to understand that the competition between short range 
magnetism and the formation of the Kondo singlet was 
an important key to understand the emerging phases in 
those compounds. For the Kondo Breakdown theory the 
reasoning is similar ; two forces compete at intermediate 
energy scales, a frustrated magnetic force leading to the 
formation of the spin liquid and and the Kondo interac- 
tion leading to the formation of the heavy Fermi liquid. In 
the phase diagram in Fig [2j To is the temperature above 
which the entropy R In 2 is released. In the Kondo Break- 
down theory, To is understood as the bandwidth of the 
f- spinons, and thus is associated with quenching of the 
entropy through the RKKY interactions, with formation 



of a spin liquid. The formation of the heavy Fermi liquid 
occurs further on the right side of the phase diagram, 
where the conventional heavy fermion metals are loca- 
ted. For convenience, it is not represented in Fig [2] It is 
interesting to notice that the two-fluids model of Ref. [49j 
comes to the same conclusion, that is in the compounds 
for which anomalous transport and thermodynamic pro- 
perties have been observed-, the formation of the spin 
liquid occurs before the formation of the heavy Fermi li- 
quid. Likewise in the vicinity of the Kondo Breakdown 
QCP, the temperature To is associated to the RKKY in- 
teractions rather than to the mean-field Kondo scale. 



THERMOPOWER IN THE VICINITY OF A QCP 
In the vicinity of the SDW QCP 

We now turn to the study of the thermopower in the 
vicinity of a SDW QCP. The quantum critical regime 
is described in Ref. [29j while some insight about the 
saturation in the zero temperature regime can be found 
in |17| . Here we summarize these two bodies of results 
and present a derivation of the thermopower in all the 
regimes around the QCP. We start with the definition of 
the thermopower as a ratio of two correlation functions 



S = 



eTLn 



(7) 



where L12 is the correlation function between the heat 
current and the electrical current, and Tn is the current- 
current correlation function defined as 

L n = linwo h Im Jo dre " T ( T -j( r ) ' -K )) - 

L u = lim^o i Im jj dre^ (T T j Q (r) • j(0)> , 

where jq is the heat current and j the electric current. 
Those two operators can be put into the following form 



(8) 



where v p = de p /dp is the velocity of the quasi-particles, 
that we consider un-renormalized by the fluctuations, and 
A(p, wjjs the spectral function. We use here the notation 
of Ref. 1 29] and define it as 



(u;/^-ep) 2 + T p 2 H 



Here Z u is the quasi-particle weight defined as 

dReYi c {kp,uj) 



z; 
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Figure 2: Tentative phase diagram of the Kondo Breakdown QCP in the presence of magnetism. The third axis 
represented on this diagram is the axis of "frustration". It can be any external parameter which competes with the 
AF long range order. When the frustration parameter is strong enough, AF disappears, revealing the Kondo 
Breakdown QCP. Within this 3D phase diagram, one observed a line of Kondo Breakdown QCPs, which are 
uncorrelated with the magnetic order. The crossing of the two critical lines of AF LRO and Kondo Breakdown is 
accidental. In the Kondo Breakdown theory, the compound YbRli2Si2 is situated at the crossings ; CeRhlns would 
be situated somewhere on the frustration axis, URu2Si2 would be deep in the heavy Fermi liquid phase (with a 
super-conducting instability at low temperatures) and CeCug-^Auj; is located at the AF QCP of itinerant 
character. This phase diagram suggests that the Kondo Breakdown QCP is a generic feature of any heavy fermion 

phase diagram; it is a universal fixed point, of non magnetic character, whose influence on transport properties 
dominates other scattering mechanisms in the quantum critical regime. Note that another phase diagram has been 
proposed [48] where the crossing of the Kondo breakdown line and the AF line has a finite width. 



and Tp(u) is the transport scattering time, which in- 
cludes both the effects of the impurities and the scatte- 
ring through the fluctuations of the bosonic mode. Ano- 
ther difference with Ref. [2j| is that r p depends on the 
position of p on the Fermi surface. We use the Mathies- 
sen's rule for adding the resistivities to get 



imp 



L(P' W ) 



dyn 



(p,<*>) . 



(9) 



To simplify the discussion, we take T^ p (p, ui) = Tq 1 as 
a constant of p and oj. The elastic scattering time tq 
encompasses for example the scattering through impurity 
centers. The effect of the fluctuations are described by 
T dyn(P^) = T h 1 in the hot regions and T^ ct (p,u;) = 
TjT 1 in the cold regions. Typically in the SDW theory 
the inelastic part of the scattering time has the following 



form 



A c . T 2 , 



(10) 



where Ah and A c are non universal constants, t^ 1 has 
the typical Fermi liquid exponent while t^ 1 has an ano- 
malous exponent due to the scattering through the soft 
quantum modes present at the QCP. Details of the eva- 
luation of Ln and L\2 can be found in the Appendix. 
The result is : 

V h V c 



L u = 



nv 2 F p* 



(11) 



where p*de = Jq + °° p 2 dp/(2ir) 2 and Vh (resp. V c ) is the 
volume of the hot (resp. cold ) regions of the Fermi sur- 
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face, satisfying Vh+V c = Vp, the total volume. For an AF 
in D = 3 where we take a spherical Fermi surface with hot 
lines at the angle O we g et Vh — sin^ A</>(T) ~ \/T 
and V c = 2 — sin0 o A(f>(T). In the case of two dimen- 
sional fluctuations in a 3-D metal, as in Ref.[29], a full 
portion of the Fermi surface is hot, even at zero tempe- 
rature. In that case Vh and V c can be taken as constants 
of the temperature. The formulae is typical of elec- 
trical transport around a SDW QCP. It can be unders- 
tood in the following way. At zero temperature, the re- 
sistivity saturates to the value Lu — •kv f PqVftq/2. At 
very low temperature for which rZh < t _1 . the correc- 
tion to the residual resistivity acquires an anomalous ex- 
ponent Ln = ttv'p Pq Vfo t// ^ Tq / 2 . Note that although this 
exponent is universal, its regime of validity can be quite 



small, since it requires that T < (r 1 ) 



-l\2/(d-2) 



A good or- 



der of magnitude for the validity of this regime is that the 
variation of the resistivity p — po (or of the conductivity) 
over which this regime is observed must be of the same 
order of magnitude as p$ itself. At even higher tempera- 
ture, the resistivity is short circuited by the conduction 
electrons, leading to a typical form Lu = %v%£nVpT c /2. 
These results are described in details in Ref. [44]. 

Let's now treat the off-diagonal correlation function 
between the heat current and the electric current. Here 
too, we have two contributions, one from the hot part of 
the Fermi surface and one from the cold part. From Eqn. 
|8]we see that L12 is odd in frequency. For the contribu- 
tion not to vanish, some asymmetry has to be introduced 
either in the summation over p via an asymmetry in the 
density of states or in the summation over uj via an asym- 
metry in the scattering times. For this purpose we make 
the phenomenological assumptions for both the impurity 
scattering time and the scattering time of the electron 
over the quantum critical modes. 7^ = t^ 1 +Lup' Ai mp 
and Tjy n — t^}(1+ta oj) where ta represents the asym- 
metric part of the scattering rate ; it has the dimension 
of a lifetime. As for Lu we find two contributions to L12 
coming from the hot and the cold regions of the Fermi 
surface. 
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hot regions of the Fermi surface while r^ 1 = 7mE£ is 
the inverse scattering time in the cold regions. From \Y1\ 
we can see that in the QC regime, the contribution from 
the hot lines is dominated by 1/Z U since this quantity 
diverges as T( d ~ z ^ z . Considering that in the quantum 
critical regime Lu saturates in Eqn. (fTTj) we get the fol- 
lowing asymptotic form in the QC regime 
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where So is the saturation value of the Seebeck coefficient 
at zero temperature. This result is quite remarkable since 
it shows that for all configurations of the hot lines, the 
correction to the thermopower divided by the tempera- 
ture tracks the variation of the Sommerfeld coefficient. 
Indeed, when the hot region has finite width, So — > 00 
and S(T) ~ 1/Z U . When the hot regions have the shape 
of a line or a point V h ~ VT and l/Z* - j 1 ( d -3)/2 so 
that the product tracks the Sommerfeld coefficient. This 
result was obtained in |29| in the case of two dimensional 
fluctuations in a 3D metal [however in this case the ther- 
mopower diverges at the QCP since the hot region has 
finite width, which is not the case when the hot region 
has the shape of a line or a point] . It is quite remarkable 
that it generalizes to all cases. 

Let's examine in more details the zero tempera- 
ture regime around the QCP. In the saturation re- 
gime, L12 is dominated by the two first terms in 
the brackets, both in the hot and the cold re- 
gions. The form of Ai mp is taken from (l7| Ai mp — 
r^/po (1 - {npoZoU) 2 ) I (1 + {npoZoU) 2 ) with U the 
scattering potential of the impurities. Within these no- 
tations we find that in the very low temperature regime 
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ImYj^ is the inverse scattering time in the 



We see that typically, the sign of S/T at low tempera- 
tures is determined by the sign of the derivative of the 
quasiparticle density of states at the Fermi level. For a 
typical heavy fermion compound, the hybridization bet- 
ween the f and c electrons lead to the following density 
of states p*(e) ~ 2p D/\if — e\ where po is the conduc- 
tion electron density of states, D is their bandwidth and 
if = gf — Sy(0) is the potential of the re- normalized f- 
levels [17j. We understand as well from Eqn. (]14[) that for 
the Yb- based compounds the S/T is negative since the 
Kondo resonance, described here by the f-level, lies be- 
low the Fermi energy (the Yb atom having 13 f-electrons, 
the shell is almost full) whereas the Ce-based compounds 
have a positive S/T since the f-level lies above the Fermi 
energy (the Ce atom has 1 f-electron, so that the shell is 
almost empty) . In the vicinity of the QCP, the residual 
I S/T I is dominated by \/Z$ and we expect it to have a 
maximum at the QCP. 

The topology of the Seebeck coefficient in the vicinity 
of a SDW QCP is summarized in Fig. [3] The most no- 
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ticeable fact about it is that \S/T\ is symmetric around 
the QCP. It diverges at the QCP in d = 2 and increases 
like T 3 / 2 and then saturates in D = 3 (see Table [TJ) . As 
a universal feature, the variation of |<S/T| with tempe- 
rature follows the one of the Sommerfeld coefficient 7. 
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Figure 3: Generic form of the Seebeck coefficient for 
the SDW theory in the case where the hot region has 
the shape of a line or a point. A saturation is present at 
the QCP, which is not the case for 2 D modes in a 3 D 

metal, as described in Ref. [29j . The left side of the 
phase diagram corresponds to the AF phase while the 
right side corresponds to the paramagnetic phase. The 
various curves correspond to different value of the bare 
mass, which describes the proximity to the QCP. At the 
QCP m = 0, in the AF phase m < while in the 
paramagnetic one m > 0. We observe a divergence of 
\S/T\ at the QCP, which is a generic feature 
independent on the presence of the hot lines. Note that 

in the AF phase, the saturation is a little bit more 
abrupt than in the paramagnetic phase. Apart from this 
small asymmetry feature, the phase diagram has a 
rough symmetric character, typical of the SDW phase 
transition. 



In the vicinity of the Kondo Breakdown 

The main difference between the SDW scenario and 
the Kondo Breakdown resides in the localization of the 
f-electrons, in the spin liquid side of the transition. In the 
Kondo Breakdown scenario, the quenching of the entropy 
is done through the formation of the spin liquid. We can 
say that the spins become entangled with one another, 
due to the presence of either the geometric frustration, 
or the frustration generated by the RKKY interactions. 
In the spin liquid phase the f-electrons are not entropic 
carriers anymore, and their contribution to the thermo- 
power is negligible around the QCP. One can see in figure 



Uthat around the Kondo Breakdown QCP, \S/T\ shows 
a pronounced asymmetry. On the right-hand side of the 
phase diagram, which corresponds to the formation of the 
heavy Fermi liquid, |<S/T| shows the same generic struc- 
ture as in the SDW case. The main response is carried by 
the conduction electrons, and the scattering through the 
QC modes is dominant. From Eqn. (|12|) . 1/Z U is given by 
the scattering through the critical bosonic modes (here 
corresponding to the condensation of the holons which 
form the heavy quasi-particle). Away for the QCP, we 
observe an increase of \S/T\. followed by a saturation at 
lower temperatures. As we come closer to the QCP, the 
value of the saturation increases, until it diverges at the 
QCP. 

A crucial difference with the SDW scenario is that the 
critical modes have a dynamical exponent z — 3 and not 
z = 2. As a consequence |5/T| diverges now in dimen- 
sion d < 3, with respectively a sub-logarithmic exponent 
j>-i/3 j n ^ _ 2 anc [ a logarithmic variation Log(T /Tq) in 
d = 3. We believe it is possible to detect experimentally 
the difference between the z = 3 regime of the Kondo 
Breakdown and the z = 2 regime of the SDW scenario, 
as will be developed in the next sections. Hence in the 
QC regime, the variation of \S/T\ with the temperature 
tracks the one of the Sommerfeld coefficient [note that 
the whole Fermi surface is hot in this easel, with 
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The most interesting observation concerns the left 
hand side of the phase diagram, where the f-electrons 
have localized. As we said, they cannot participate any- 
more in carrying the entropy, which leads to a dramatic 
discrepancy from the SDW phase diagram. In this part 
of the phase diagram, if the AF order is present, a si- 
gnature milder than for the SDW should be observed in 
\S/T\. The change of \S/T\ as we pass through the AF 
transition is typical of the proportion of itinerant versus 
localized character of the magnetic order. If the magnetic 
order comes mainly from localized electrons, the jump in 
the thermopower coefficient at the transition should be 
mild, due only to the indirect opening of a gap in the 
critical modes, as a consequence of the formation of the 
order. On the other hand, it the magnetism is due to the 
formation of wave from itinerant electrons, we can expect 
that the response in the thermopower will be significant. 

The thermopower is very sensitive to the scale E* of 
the Kondo Breakdown. The scale E* is the energy be- 
low which the mismatch between the conduction elec- 
trons and the spinons Fermi surface becomes noticeable. 
Above this energy scale the Seebeck coefficient follows 
the QC regime ; it is dominated by the scattering of the 
conduction electrons through the QC modes (here the 
holons). Below E* , the propagator of the QC modes is 



gapped, it means that the effective hybridization between 
the conduction electrons and the spinons vanishes, and 
the Seebeck coefficient is dominated by the residual Fermi 
liquid contribution. This results in a dramatic drop of 
\S/T\ below E* , since when the temperature is decrea- 
sed, the conduction electron's scattering changes brutally 
from the QC to the Fermi liquid regime. The two regimes 
are physically disconnected from each other. This brutal 
drop at E* is similar to the brutal decrease expected in 
the T = limit when we cross the QCP from the heavy 
Fermi liquid side towards the spin liquid side. As seen in 
Figure[3]thc reconfiguration of the Fermi surface is revea- 
led by a brutal drop of |SyT| when going from the Fermi 
liquid to the spin liquid phase. Whether the Seebeck co- 
efficient changes sign or not will depend on the details of 
the conduction scattering in the spin liquid phase and will 
vary form compound to compound. A strong signature of 
the Fermi surface reconfiguration should be observed in 
the localized phase. 

We keep in mind here that a weakness of the theoretical 
treatment of the Kondo breakdown model reside in the 
fact that the localized degrees of freedom are badly taken 
into account at the present stage of the theory. However 
the Seebeck coefficient is precisely less sensitive to those 
degrees of freedom, since they don't carry entropy. That's 
why it is maybe a decisive test for revealing the scale E* 
and thus differentiate between the SDW and the Kondo 
Breakdown QCPs. 

There is a potential issue of whether the scale E* can 
be masked by the occurence of AF order. To be more 
precise, insensitivity of the Seebeck coefficient on the AF 
order is observed near the QCP, if the energy scale E* 
is larger than the Neel temperature Tjv- Away from the 
QCP in the AF side, in the case where Tm > E* , the si- 
tuation is somewhat complicated. The whole "Fermi sur- 
face" of spinons can be gapped below the Neel tempera- 
ture, and holon excitations will become already suppres- 
sed even above E* . In that case, the abrupt drop in the 
Seebeck coefficient might occur from Tjv instead of E* . 
This feature is also completely different from the SDW 
scenario. If cold regions are still present at the QCP, there 
will the measurable signature in the Seebeck coefficient 
at the Neel temperature should be much milder than the 
abrupt dropping at E* which occurs in the Kondo break- 
down model. 



A SMALL SURVEY OF FOUR COMPOUNDS 
YbRh 2 Si 2 

The most recent results of the thermopower in 
YbRli2Si2, driven to the QCP via the application of a 
small magnetic field, comes form Ref.[13|. For an applied 
magnetic field B < 65mT , a negative Seebeck coefficient 
S < is found, in good agreement with other measure- 
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FIGURE 4: The form of the Seebeck coefficient divided 
by the temperature around the Kondo Breakdown QCP. 
The diagram shows a pronounced asymmetry between 
the heavy Fermi liquid phase and the spin liquid phase 
[56|. \S/T\ shows a brutal drop at T — when going 
form the Fermi liquid phase to the localized phase, 
which reveals the re-configuration of the Fermi surface 
at the QCP. The scale E* signals a brutal change from 

the QC to the spin liquid regime. S/T seems to be a 
very good experimental probe to detect the multi-scale 
character of the Kondo Breakdown QCP. 



ments for Yb compounds (Hi - 54 1 . A logarithmic increase 
-S/T ~ - log(T/f ) with f = 3K is observed in the QC 
regime, which is defined for this compound as the regime 
for which B — Q5mT and T < 25K. This logarithmic law 
is observed above a temperature T max = O.IK. Below 
T m ax , —S/T drops abruptly to reach a very low value. 
In Ref. [28] a change of sign is associated with this abrupt 
drop of the Seebeck coefficient, and it is argued that at 
B = 0, in the left hand side of the phase diagram, the 
Seebeck acquired a positive value. It is to be noticed that 
at the QCP, above T max the variation of —S/T follows 
the variation of the Sommerfeld coefficient. This behavior 
changes below T max with the sudden decrease of —S/T 
whereas the Sommerfeld coefficient shows an upturn as 
the temperature is lowered. Another very anomalous pro- 
perty is that no sign of the magnetic phase transition is 
observed in —S/T: the only temperature scale observed 
in this part of the phase diagram being T max , the scale 
of the abrupt decrease. 

This body of results can be simply interpreted with the 
Kondo Breakdown theory. The abrupt change in —S/T 
is attributed to the reconstruction of the Fermi surface 
around the Kondo breakdown QCP. The logarithmic in- 
crease is naturally interpreted with a z = 3 QC regime, 
in d = 3, which is precisely the prediction of the Kondo 
breakdown theory. The evidence for a z = 3 QC regime in 
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this compound, is supported as well from the logarithmic 
variation of the Sommerfeld coefficient with temperature 
and the variation of the Griineysen ratio like T~ 2 / 3 [42l|. 
The phase for B < can be interpreted as the phase 
where the f-electrons localize, and T max is within the KB 
theory the scale E* below which the conduction elec- 
trons become insensitive to the scattering through the 
QC modes. The discrepancy below E* between the va- 
riation of —S/T and the Sommerfeld coefficient 7 seems 
to indicate that the upturn in 7 is due to the presence of 
localized moments, which don't participate in the trans- 
port of entropy. 

In Figure [2] we have put YbRli2Si2 exactly at the cros- 
sing point between the AF line and the Kondo Break- 
down line of QCPs. Within the interpretation of our 
theory, it is just an accidental fact, but this interpreta- 
tion is under debate within the community. It would be 
very interesting to get an experimental insight on what 
happens when one goes a bit away from this intersection. 
Recently, YbRli2 Si2 has been doped with a few % of Ir 
and Co. The Ir doping pushes the compound outside the 
AF phase while the Co doping pushes it inside the AF 
phase. It would be of the greatest interest to measure 
the thermopower in the case of Ir and Co doping. Within 
the Kondo breakdown theory we expect that the features 
of the pure compound will be reproduces with no major 
changes as soon as E* > Tjy. The presence or absence 
of magnetic order should not affect in a major way the 
location of the scale E* ( called T max in Rcf[28]). For do- 
ping with Co, we will probably be in the situation where 
E* > Xjv, in which case we can expect some changes to 
start at Tjv rather than E* but those will be of small ma- 
gnitude compared to the drop at E* . The abrupt drop 
of the thermopower at the QCP can indicate in a precise 
manner the location of the Kondo breakdown QCP, and 
as such the measurement of thermoelectric effects is de- 
cisive in corroborating or invalidating the results of Ref. 
[HI . Of particular interest is the confirmation of a QCP 
under the AF dome in the Co-doped YbRb.2Si2. Lastly, 
some measurements under hydrostatic pressure are ne- 
cessary to validate the whole picture [5 7]. 

In this compound the application of pressure or do- 
ping has important consequences on the structure of the 
magnetism. It mainly affects the amount of frustration 
and the dimensionality of the magnetic order. The ther- 
mopower however, occurs to be mostly insensitive to the 
details of the magnetic order, as soon as it is of localized 
character. That's what makes it such an attractive expe- 
rimental probe to test the Kondo breakdown scenario. 

Ce(Mn)In 5 

In this series of compounds, the only study of the ther- 
mopower close to a QCP concerns CeCoIns [26]. The 
compound is superconducting at T = 0, and driven to a 



QCP around H c2 ~ 5.5T. Around this field driven QCP 
the thermoelectric properties have been thoroughly in- 
vestigated in Ref. [26l|. The results certainly don't show 
any strong re-structuration of the Fermi surface at the 
QCP. When the applied field crosses the QCP, a small 
increase of S/T is observed. The Seebeck coefficient is po- 
sitive in the whole phase diagram. These results suggest 
that the QCP is not a Kondo Breakdown, but maybe 
in the universality class of an AF SDW, or in a third 
universality class, that could appear in the presence of 
a strong magnetic field. In particular, the fact that the 
ratio q = (SNAe)/(-fT) departs from unity at the QCP, 
which is well reproduced by the SDW scenario [l7j . 

In Figure [2] we have put CeCoIn5 on the right side of 
the Kondo Breakdown QCP (the SC phase has not been 
represented here). It is possible that under a magnetic 
field, CeCoIns is driven towards a QCP associated with 
short range AFM. 

For this series of compound, the best change to find the 
Kondo breakdown QCP is around CeRhliis. This com- 
pound is an AFM at low temperatures. With the appli- 
cation of pressure, it is driven towards a phase transition 
around l.75GPa where the Fermi surface re-configures 
[58| . It our belief that this QCP is associated with the 
Kondo Breakdown, with a quasi-two dimensional nature 
of the QC fluctuations[59]. If it is the right hypothe- 
sis, the measurement of the thermopower under pressure, 
around the point where the Fermi surface re-configures 
should show a dramatic change, with S/T dropping off 
from the heavy Fermi liquid at high pressure to the local 
f-electron phase at low pressure. An interesting point to 
investigate here is whether the magnetic order at zero 
pressure is of itinerant or localized character or maybe 
both at the same time(6p|. If the magnetic moments are 
fully localized character, no strong signature of Tjv shall 
be observed in S/T, whereas, the scale E* precursor of 
the reconfiguration of the Fermi surface, shall be obser- 
ved instead. On the other hand, if the magnetic order is 
of itinerant character, a signature of Tjv comparable to 
the one observed in the specific heat measurement shall 
be observed. It is possible to apply magnetic field as well, 
leading to a line of QCP with the magnetic field. If it is 
possible, it would be very interesting to measure the ther- 
mopower close to this cline of QCP. When following the 
line of QCPs, a cross over towards SDW type scenario 
might be observed, similar to the one found in CeCoIns. 

In this compound as well, the study of thermoelectric 
properties would be a decisive experiment in order to 
elucidate the nature of the QCP under pressure. 

CeCus-zAiia; 

This compound is one of the first where the presence 
of a QCP was detected[3|- When this compound is do- 
ped with 0.1% of Au one reaches an AF QCP. In the QC 
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regime the specific heat was shown to vary logarithmi- 
cally with temperature, while the resistivity is linear in 
T. Neutron scattering experiments @ have revealed that 
the dynamic spin susceptibility has a pronounced two di- 
mensional character, and shows anomalous exponents for 
a wide range of q-vectors in the Brillouin zone. 

Two theories are in competition for this compound. In 
Ref.0| it was argued that the QCP is very anisotropic, 
and its nature is one of a two dimensional SDW in a 
three dimensional metal. This theory reproduces the li- 
near corrections to the resistivity, as shown in Table HI 
Since the chemical doping with Au introduces a substan- 
tial amount of disorder, it is conceivable that the linear 
resistivity observed in this compound is due to a wide 
tail correction to the residual resistivity, in the d = 2 
SDW scenario. In d = 2 the SDW scenario is believed not 
to produce any anomalous exponent when the electrons 
have been integrated out of the partition function (HH- 
It has been argued, however that it is no so when the 
electrons are treated self- consistently with the quantum 
critical modes. In that case, anomalous exponents have 
been predicted for the staggered dynamical spin suscep- 
tibility, in the QC regime [62| . The observed linear varia- 
tion of the transition temperature Tjv with the doping x 
corresponds as well to the two dimensional character of 
the QC fluctuations within the SDW theory. 

Another theory has been proposed to explain the ano- 
malous properties, called the locally quantum critical 
scenario [47J. This theory assumes that the bosonic modes 
have two dimensional character, and then using extended 
Dynamical Mean Field Theory (e-DMFT) it is argued in 
this work that a local mode emerges at the QCP, leading 
to a reconfiguration of the Fermi surface. 

As a result the right theory for this compound is still 
very mysterious ; is it a SDW with two dimensional cha- 
racter, or a more unconventional locally quantum critical 
scenario ? 

The only study of thermoelectric effects is a short 
note where it is shown that the thermopower diverges 
at the QCP [24| . It is not clear what is the exponent 
of this divergence. It would be very interesting to mea- 
sure the thermopower in the AF magnetic phase of this 
compound. If the Fermi surface is re-configured around 
the QCP, an abrupt change of S/T is predicted. On the 
other hand, if the two dimensional SDW scenario is the 
right answer, one should observed in S/T a consequent 
signature of the Neel temperature or the order of the one 
observed in the Sommerfeld coefficient 7. Moreover, the 
phase diagram in that case will be show broadly symme- 
tric features between the ordered phase and the parama- 
gnetic phase, with no abrupt change of S/T at the QCP. 
In Figure [2] we have placed CeCu6-xAu x at the proxi- 
mity of the AF order, indicating that in our opinion, the 
most likely theory to apply here is the one of Refe 0,123 . 
However, itis only an opinion, and we don't have enough 
substantial scientific arguments to corroborate it at the 



moment. Here again, the measurement of the thermoe- 
lectric effects appears to be a decisive experiment. 



URu 2 Si 2 

This compound constitutes one of the most enduring 
mysteries in the field of strongly correlated electrons. A 
very well defined phase transition carrying more than 
40% of the free electron's entropy occurs below To = 
17 K. Despite almost twenty five years of experimental 
investigations, the mystery concerning the nature of this 
"hidden order" remains unsolved [63] . Several theoretical 
proposals have been made, including some exotic short 
range anti-ferromagnetism [6~i| , a Lifshitz transition [65| , 
a charge density wave scenario [75| and a scenario where 
the hidden order has a localized character 166 
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experiments [68j have revived the bad structure studies. 
A debate exists on the nature of the two f-electrons in 
the U-f 2 atom, whether the localized picture is correct 
[69l | or whether the itinerant one [70] is the correct pic- 
ture. Very interesting experiments under pressure show 
a long range AF (LRAF) order occurs at the pressure 
of 0.5GPa, and the transition form the hidden order to 
the LRAF order is of the first order 7l|. Neutron scat- 
tering experiments show two types of excitations, one at 
Qaf = (0,0,1), which becomes static and long range 
under pressure in the AF phase [72^. and another excita- 
tion at Qo = (1, 0, 0), which is present only in the hidden 
order phase[73]. The superconducting phase is as well of 
very unconventional nature [74] The study of thermoelec- 
tric effects in URu2Si2 is complicated by the fact that it 
is a compensated metal. Below To it has been established 
that the number of carriers drops considerably, leading 
to the physics of very low density of electrons [76l478j . 

With all these observations in mind, it might look sur- 
prising to test URu2Si2 as a potential candidate for the 
proximity to a QCP. A recent study of thermoelectric ef- 
fects on this compounds might change this perspective. 
In the paper Ref [27J , a thorough study of both resistivity 
and thermopower has been conducted under magnetic 
field. This study first confirmed the strong anisotropy of 
this compound. The anisotropic nature of this compound 
was known for a long time, with the observation of aniso- 
tropy factor of 3 to 5 in the resistivity l79f, the magnetic 
susceptibility [80j and the critical field [8 1|. The de Haas 
van Alphen study for this compound captures only a mild 
anisotropy in the three Fermi surfaces observed. Very in- 
terestingly Ref(27| reveals a significant anisotropy in the 
inelastic scattering of the normal phase. When a magne- 
tic field of 12T is applied, anomalous scattering is obser- 
ved in the electrical resistivity of the basal plane, which 
is then linear in temperature, while the c-axis resistivity 
remains Fermi liquid like down to very low temperatures. 
The Seebeck coefficient divided by temperature shows as 
well a very anomalous behavior. In the basal plane it 
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departs from the constant value predicted by the Fermi 
liquid theory to finally change sign for fields larger than 
12T, at temperatures lower than T c i iange — 0.8K. In the 
c-axis, the signal remains Fermi liquid like, with a well 
defined saturation for all fields considered. This anisotro- 
pic situation is very reminiscent of the case of CeCoIiis, 
for which as well an inelastic transport time has been 
revealed[82j . 

This body of observations motivates us to suggest 
that URu2Si2 might be in the proximity of an anisotro- 
pic Kondo Breakdown QCP. The anisotropic scattering 
is cut-off below the scale E* characterizing the Kondo 
Breakdown theory. In the present case E* would be ani- 
sotropic, with a small value of the order of E* ah ~ 0.8K 
in the basal plane and with a much bigger value of E* 
in the c-axis direction. A detailed exposition of this pro- 
posal will be published elsewhere[83]. At the moment we 
would like to suggest that it would be extremely inter- 
esting to test this idea by exploring the thermopower on 
the whole pressure phase diagram. If the AF order is of 
localized nature, no significant entropy transport should 
be associated with the occurrence of the AF phase, and 
the signature in the Seebeck coefficient should be minor 
near the QCP. On the other hand, if the AF order is 
of itinerant character, a strong signature in S/T is to 
be observed. Moreover, if the compound is sitting at the 
proximity of the Kondo Breakdown QCP, one expects to 
see some evidence of the scale E* at other pressures in 
the phase diagram, and especially in the LRAF phase. In 
the Figure [2] we have placed URu2Si2 in the proximity, 
but still a little bit away from the Kondo Breakdown 
QCP. For this compound again, thermoelectric studies 
under pressure could be decisive to unveil the mysterious 
nature of the hidden order. 



CONCLUSION 

The aim of this discussion paper is principally to en- 
courage new experiments using the thermopower as a tes- 
ting probe for discriminating the nature of QCPs in heavy 
fermions. It turns out that the two main classes of QCPs 
in heavy fermions have very different signatures in terms 
of the Seebeck coefficient. The SDW scenario has See- 
beck coefficient with a good degree of symmetry around 
the QCP, between the ordered phase and the parama- 
gnetic phase. On the other hand, for the Kondo Break- 
down QCP, the Seebeck coefficient shows a pronounced 
asymmetry around the phase diagram, dropping out in 
the Kondo broken phase, since the f-electrons are not 
available anymore to carry the entropy and the quantum 
critical scattering of the conduction electrons is gapped 
below an energy scale E* . 

The Seebeck coefficient can be used as a very sensitive 
probe to detect whether the magnetism is of localized or 
itinerant character. In the case of itinerant magnetism, 



S/T is qualitatively tracking the variation of heat at the 
magnetic transition. For magnetism emerging from loca- 
lized moments, the specific heat is expected to be one 
order of magnitude more sensitive to the phase transi- 
tion than S/T, since in that case the localized f-electrons 
don't participate to the heat transport whereas their en- 
tropy is locally quenched by the apparition of the order. 

Lastly, the temperature dependence of S/T in the QC 
regime is tracking the variation in temperature of the 
Sommerfeld coefficient, which enables to make the dis- 
tinction between different classes of QCP, with dynamical 
exponent z = 2 or z = 3. 

It is our belief that new experiments within this tech- 
nique, especially under pressure, can shed light on the 
nature of the various QCPs of heavy fermion compounds. 
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Appendix : derivation of eqn. (|13p 

In this Appendix we derive the equations leading to 
CD}, dH and CEO). 



The current-current correlation function 

We start with 

in=Ept$r 0O du,(-^)^(p,a,) 

with A(p,u) = - 

and t^(u) := t^ 1 + r^(u) . 

The value of Tq 1 and T^ n are given in the text in 
Eqns.© and (fTO)) . The scattering time r is conside- 
red here as valid respectively in the "hot" and "col- 
d" region, and the subscript has been omitted. Using 
^ p = J D D p(e)de and noticing that the wave vector v p 
is pinned at the Fermi surface, we get 
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over e we get 



Ln = v 2 F p f™ x (-15) cfe 
changing variables for e — >■ y we get 
in = ^Po ^ (-$£) r(w) dy 



and remembering that depending on the region 

in the Fermi surface, Tdyn = i~h or Td yn = t c , we get 



P* Vh 



PoVc 



2 Vr^+r,; 1 ro^+rc- 1 



(A.16) 

with in spherical coordinates p*(e) — p^Wdi » = 
sin </>o A</>, V'c = 2 — sin (froAcj) and A</> ~ V^T the width of 
the hot regions. 



The heat-current correlation function 

The heat-current correlation function is given by 

L 12 =r P ^r^(-i£)^ 2 (P,") 



with A(p, oS) — 



(u/Z u -e p ) 2 + T-*(u) 



r X ( w ) = T imp + W« 



and ^( W ) = T Vc ( X + T A ^) ■ 



The definition of r h A is given in the text in Eqn. ([TO]) . 
The heat-current correlation function is evaluated with 
the following steps. Transforming £h/p m t° an integral 



L12 = v%V h/c / de (pq + epoO 



/ , / 5/ 
x / — — — I a; 



permuting the integrations and changing variables 
e — » £ + w /Z^ , we get 



L 12 =vlV h/c j^du:(- d J- 
x / &(pg + ( £ + c,/Z„)pg') 2 

then changing variables y :— st(uj) 

L12 = v%V h / c J duj (- ^\ ut(uj) 

x J™ dy (p* + (y t-V) + u/Z u ) Po') (^W) 

the term linear in y vanishes and the contribution 
in front of pg reads 

using the definition of t(uS) we have 

L 12 = —^j- V h/c Po(~ J J- 



the contribution in L\2 in front of pg' reads 



-12 



which gives 



L 



12 



6 Jr4 „*Vt 2 ^^J__ 1 

2 T W 



finally 

-^12 = ^12 + L 



(A.17) 
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